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1„  Introduction  and  Sn— 

In  this  paper  wo  describe  oognential  Ufa  toot  procedure*.  As 
in  a recant  paper  dented  to  non  oaqoitlal  oathods  we  consider  the 

special  caso  in  which  the  underlying  distribution  of  the  length  of  life  ia 
given  by  the  eapeoantlal  density 


(1) 


f(x,0) 


x >0 


where  the  unknown  pa  ran  *•+-—?  6 > 0 can  bo  thought  of  physically  as  the  mean 
life.  Our  primary  aim  ia  to  teat  the  simple  hypothesis  9 * 0^  against 
the  nimple  alternative  E^*  9 - (9#  > S^)  with  type  I end  II  errors  equal 

to  predesigned  values  a and  g respectively.  The  test  is  carried  out  by 
drawing  n items  at  random  fres  the  population  and  placing  them  all  on  a life 
test.  Vs  consider  both  the  replacement  case,  in  which  failed  ites*  are 
immediately  replaced  by  new  itens,  and  the  non-repiscemant  case. 

An  interesting  feature  of  the  tests  is  that  they  con  be  terminstad 
cither  at  failure  tises  with  rejection  of  K#  nr  at  any  time  between 
failures  with  acceptance  of  HQ.  Since  abnormally  long  ir terrain  between 
failures  furnish  "in  formation"  in  favor  of  Hc>  and  abnormally  abort  intervals 
furnish  "inforaationn  in  favor  of  the  above  features  are  not  only 
reasonable  but  actually  uesi ruble 0 These  features  were  pointed  out  in 
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£~ U,J  - Similar  prohlnsa  invalving  a coKilsaous  tl—  para— tar  have 
recently  appeared  In  C 'D % 

In  this  paper  we  obtain  likelihood  ratio  testa,  giva  fornuit*  for 
tha  O.C.  curve,  fbr  the  expected  mater  of  failure*  E^(r)  and  tha  axpactad 
waiting  time  EQ(t)  before  a decision  it  reached.  In  the  replacement  case 
whore  the  maker  of  its—  on  tost  throng  out  the  asperlmcit  is  the  same, 
namely  n,  it  is  shown  that  E (t)  " - EQ(r)*  A table  giving  values  of  L(9), 

9 n I, 

Ea(r)  and  Efl(t)  for  certain  choices  of  — , a,  and  0 is  given  for  the 
replace—  nt  ea 


2,  Basic  formulae 

Kald ’s  work  an  sequential  analysis  can  be  used  virtually 

without  modification  in  a situation  where  decision*  .irs  made  continuously. 
In  fact,  in  a truly  continuous  situation,  Wald’ , formula#  become  sxact 
since  there  Is  then  no  excess  over  the  boundary.  It.  will  beccas*  clear  as 
we  proceed,  that  in  the  problem  at  hand,  the  situation  can  be  termed  send- 
continuous  (not  to  be  confused  with  the  concept  of  the  same  nams  in  real 
variable  theory)  since  there  is  no  excess  over  the  boundary  used  for 
accepting  Hq,  but  there  will  in  general  be  so—  excess  over  the  boundary 
used  in  accepting 

Let  us  assume  that  the  ur<d erlying  pad.f„  is  (l).  n items  are 
drawn  at  random  Trm.  (1^  and  placed  on  life  test.  We  wish  to  test 
H^:  9 - 9C  against  H^s  3 « with  type  I error  » a and  type  II  * 0. 
Information  is  available  continuously  and  it  is  readily  verified  that  a 
continuous  analogue  of  the  sequential  probability  ratio  test  of  Wald  can  be 
used.  The  decision  as  time  unfolds  depends  on 


(2) 
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idiero  B and  A an  constants,  on  a and  p,  such  that  B < 1 < A. 

Th«  decision  to  csntinus  sjperlasststion  ia  Bade  as  long  as  the  inequality 

(2)  holds.  If,  at  the  tims  the  experiment  is  stopped,  the  first  inequality 

in  (2)  is  violated  we  accept  H#;  if  the  second  inequality  is  violated,  we 

fi 

accept  ft.  As  in  Wald's  case  the  test  obtained  by  setting  B * — - — , 

(1)  1 “ u 

A *>  ie  a satisfactory  solution  of  the  problem'  1 from  a practical 

Cl 

point  of  view. 

V(t)  in  (2)  is  a statistic  which  can  be  interpreted  as  the  total 
life  observed  up  to  time  t.  Is  the  replacement  case 

(3)  V(t)  - ni0 


(l)  It  lia3  een  pointed  out  by  llaid  that  in  order  to  have  a teat  of  exactly 
strength  (a,p),  A and  B in  (2)  should  be  replaced  by  A*  and  B*,  where 


In  the  present  case,  due  to  the  fact  that 


A*  < i-u  and  B*  > JL 

“ a — l-o 

information  is  available  continuously  in  tins,  wa  know  that 

e-  B 

B - B • - — , aince  acceptance  of  H involves  no  excess  over  the 
1-a  ° 

boundary.  However,  acceptance  of  H.  does  in  general  entail  s positive 

**  # 

excess  over  the  boundary*  ?nd  all  we  can  say  initially  about  A is  that 


it  should  lie  between  A9^/S#  and  A0  Thus  using  A “ 3na+e*d  of  A" 


1~B 

- » ‘ 1 nan 

a ' ~~T~ 

is  an  approximation.  The  approximate  test  based  on  u«Ing  A and  B is 
completely  suitable  since  its  strength  (a',p’)  is  such  that 


<j,  &<8*  < 


g 

1-a 


i*  ♦ 0 s < a + p /*7,  PPo  <+5 -q_/.  Sine® 


a and  £ are  generally  small  (<  .10  say),  a procedure  based  on  using  A 
and  B provided  essentially  the  same  protection  agr.inst  errors  of  the 
first  and  second  kind  as  does  a test  based  on  A*  'Old  B*.  While  cue 
can,  using  formula#  ir.  /”" 2_7,  coSiput*  A*  ir.  such  a w ay  as  to  give 
exactly  strength  (a,0),  there  seems  to  be  little  reason  in  most 
practical  problems  for  expending  the  time  and  eifert  involved. 


1 


(2' 

In  the  non-re  pi  ee — ent  case  * V(t)  is  given  tog r 


U:)  7(t) 


r r 

XT  ♦ (n-r)(t-xr)  a ♦ (n-rXt-x.), 


where  x(  denote*  the  tine  of  the  ith  failure  (x_  * 0)« 

It  is  convenient , if  one  wishes  to  graph  the  data  continuously 
in  time,  to  -rite  (2)  in  the  fern 


(5) 


* rs  < V(t)  < h#  ♦ re 


where  teQ,  h^,  and  s are  positive  cores  tint*  given  ty 


(6) 


- log  B 

J,-  i 

h *• 


"1“ 


leg  A 


1 1 


log 


®o 


1 1 


Further  it  can  be  sbowu  as  in  £"lf  pp.  i*8-50_7  that  the  0;C5  curve,  i0ec 

the  probatoijJLty  of  accepting  H when  9 is  the  true  parameter  value,  it 

o 

approximately  given  by  a pair  of  parametric  equations 

(*o\» 


(7) 


L(«) 


U 

a — x 


a - 


t T-  I - 1 
< CS  / 

' 


U v 


by  letting  the  parameter  h run  throu#i  all  real  valuee0 


(2)  In  the  non-replacement  case  it  may  happen  that  no  decision  has  been 
reach*?’  by  the  tin*  t ■ x , when  all  n i tarns  have  failed.  This  Kill 
then  re  quire  that  we  either  put  mere  items  on  test  and  wait  until  (2) 
is  vio  ited  or  else  have  a rule  which  will  tell  us  how  to  terrains!* 
the  experiment  and  with  what  decision  at  t • Fortunately  r»  is  of>n 

at  our  disposal  and  sc  can  be  chosen  sufficiently  large  30  tist  the 
probability  of  reaching  no  decision  by  time  Xj^  is  negligible,.  For  1 
enough  n,  it  really  makes  very  little  difference  how  we  truncate  experi- 
oentrtion.  We  could.  , adopt  the  rule  that  is  accepted  if  (?) 
ie  satisfied  for  all  t < x^. 


t 


4tin> 


-y 


The  valus  of  L(9)  at  the  Uve  points  9 ■ 0,  9^,  a,  0O  arsd  Ov 
enahles  one  to  stench  the  entire  curve.  Thee*  value*  are  respectively  0, 

B,  log  A/(log  A - log  B),  l-«,  and  1. 

We  now  give,  in  terms  of  L(0),  a formula  for  E0(r),  the  expected 
number  of  observations  required  to  reach  a decision  when  0 is  the  true 
parameter  value.  Since  the  logarithm  of  the  middle  expression  In  (4)  is 
sitter  log  B or  log  A at  the  time  experimentation  stops,  ws  have,  neglecting 
only  the  exceee  over  log  A, 

(8)  Eg(r)  log  g~  - E9(V(t))  E0(V(t))^L(9)  log  B ♦ [l  - L(9)J  log  A 

It  i*  prcrad  in  the  next  section  that 


(9) 


Ea(V(t»  - © Ka(r). 


Hence  we  have  from  (9)  and  (10) 


L(8)  lag  B fl  - L(9)]  log  A hI  - L(9)(ho  ♦ hx) 

1)  - 


'10>  F,(r)w/  w log  A log  B hh, 

l rnr  “ '■*  ■ . 


for  9 - es 


e - 


If  we  let  k “ Qq/9^,  K9(r)  be  com  particularly  simple  vhcc 
, e,  or  9#.  The  result  is 

0-1)  E^(rW0  log  B ♦ (1  - «)  log  h}j R eg  k ♦ (k  - i)j 

2 

Ew(r)^'~  log  A log  B/(log  k) 

E«  (r)<^[^(l  - a)  leg  3 ♦ a log  Aj  /&xg  k - (k  — l)j  - 


for  9^3 
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In  table  1 w*  give  Kg(r)  for  the  five  values  9-0,  9^  e,  0^,  oe  , for 
four  values  of  k(3/2,  2,  5/2,  3),  and  the  four  number  pairs  (a,0)  /hich  can 
be  made  with  the  maebers  .01  and  ,05» 

3.  A baric  Id  to  tit  y 

In  this  section  (9)  is  derived.  While  this  result  car.  be  obtained 
as  a consequence  of  a theorem  of  Dcob  on  continuous  parameter  martingale* 

.-'"l,  p.  376_7,  it  seems  desirable  to  give  & simpler  proof.  We  shall  consider 
ths  replacement  case,  where  V(t)  - nt,  although  the  proof  can  be  trivially 
modified  so  as  to  hold  in  non-replacement  and  truncated  situations,, 

In  the  replacmsent  case  (9)  becomes 

(12)  Ke(t)  - SQ(r)  | . 

Thus  we  are  relating  expected  waiting  time  to  reach  a decision  to  the 
expected  nui,  ~er  of  failureSo 

To  prsrr  (12)  vs  introduce  a "large,"  integer  N and  let 
denote  the  time  of  the  Nth  failure.  Lei.  t denote  the  (first)  time  at 
which  the  inequality  (2)  is  viol* ted  or  whichever  comes  sooner.  Then 
ve  can  write 

(li)  *,-<■*  <Vi  <Va  - W * •••  * (*n  * ’Vi5- 

We  either  accent  before  sjj,  in  which  case  t - Xj,,  or  accept  before 
Xjj,  in  vinich  case  3^  < t < x^, , pr  take  no  action  before  Xjj,  in  thi  ch 
case  t - and  r » Nc  Sine*  w is  fixed  in  advance 

(14)  EvXj,)  - N | « 

Further  it  is  easily  verified  that  the  (N-r)  random  variables 
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(x_^  - t),  (3^2  ” • ••#  (*h  “ *5-1)  ar*  u*i*P«*i*»ntly  and 

identic  ally  distributed  with  the  exponential  dena  ity  ~ e"553^',  x > 0o 
Hence  if  we  take  the  expec  titlOA  of  loth  sides  of  (14) » first  bolding  r 
fixed  and  Lfarn  taking  the  expectation  with  respect  to  r,  we  obtain 

(15)  K | - EqUjIO  ♦ f«  - %(r;»)l  2 

n *•  v j a 

•r 

(16)  Ee(t?H)  - S^r;*)  I - 

Formula  (16)  holds  for  ell  N.  Far  K— ^ the  probability  of  coming  to  a 

decision  before  x^  tends  to  unity,  Moreover  as  N — b»o 

(17)  E^rjH)  f EQ(r)  and  Ee(wN)  f EgCt), 

vAiere  E0(r)  and  E0(t)  are  respectively  the  expectsd  number  of  failure!  and 
expected  waiting  time  to  reach  a decision  if  N ■ . Thus  it  follov.-et 

letting  E — >33,  that  (16)  becoaee 

(12)  EQ(t)  » Eg(r)  | . 

The  non-replacement  caas  can  be  treated  in  exactly  the  seas  *ay0 
This  is  because 

(18)  V(a^)  « V(t)e(n-r)(xr^1-t)e(n-r-l)(xw>2-xJ^1)  + ♦ (*n~:^x)o 

Ae  before,  t * x^.,  if  % is  accepted;  x^.  < t < H*  i»  accepted; 

ana  t « 3^,  if  no  decision  is  reached  by  the  time  all  n itans  have  fail-c,. 
The  last  (rv-r)  component*  on  the  right-hand  side  of  (18)  are  mutually 
independent  randan  variables,  each  distributed  with  the  ped»ft.  (!)<•  Thus 
it  irollwa  as  in  the  r«c  placement  case  that 


(19) 


S0(7(t)jn)  - O Te(r;n); 
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As  n increases,  E^(rjn)  f Eg(r),  ths  expected  number  of  failures  in  reaching 
a decision  5xi.  the  repi&csmsnt  csss.  Thus  no  ostia  r how  vs  d sc  ids  to  tsminato 
sxpsrlosntat.isn  ®0(4 * * 7(t);n)  can  he  replaced  by  2Q{V(t))  » 9 E&(r),  n is 

largs.  In  p;  >ctiee,  for  "large"  n ons  could  taka  n > 3 ngx  EQ(r)0 
Rsaark;  It  should  bs  noted  that  while  vs  can  relate  expected  waiting  tins 
to  the  sxpsctsd  nuobsr  of  failurss  in  ths  replacement  case  by  (12),  (19) 
relates  expected  total  Ilfs  (not  waiting  tins)  to  the  sxpsctsd  men  bar  of 
failurss  in  ths  non-re  pieces*  nt  cass.  Actually  ons  has  to  know  ths 
probability  distribution  of  r in  order  to  compute  Kg(t,)  exactly  in  ths  non- 
replacement  cass.  It  can  bs  shorn,  in  ths  non- replacement  case  that  ths 

fftwwil  ta  ) i;  Ky 

(30)  E0(t)  - II  Mr  - k|e)  Ee(xk>n), 

t- 

R 

where  E0(X,  ) « 9 > l/n-i+1.  In  ths  replacement  cass  one  has 

(analegous  to  (20))  ths  formula 

(20*)  EQ(t)  -IT  Pr(r  - k|e) 

au— X 

tiler*  n is  the  sample  sizs  maintainsd  throughout  the  experiment  end 

n^  " If®/0*  Thus,  in  the  replacement  ca*s  (20‘ ) clsarly  becomes  (l?)„ 

(20)  is  valid  for  a 11^  life  test  procedures  which  involve  non -r eplacsmen t = 

(L) 

Similarly  (3C’)  holds  for  allv  life  tsst  proesdurss,  where  items  which 
fail  are  replaced, 

(4)  i^So,  for  truncated  or  untruncated,  sequential,  or  any  similar 

procedure.  Of  course  ths  probability  distribution  of  r dome  d»p*nc 

on  the  procedure  which  is  followed.  In  sxplici  t formulas  for 

Pr(r  * k}9)  are  workod  cut  for  tires  procedures. 


i'°bXe®  ; Kind  a sequential  rtplic—nt  procedure  for  beating 

’o’  9*  7500  hour*  against  H]_:  0^  • 2500  hours  with  a * 0 - c05.  The 

; on  star,  number  of  items  under  tart  i»  n • 100. 

Solution  . In  this  case  (2)  beoomer 

. r -7{t)/3?50 
53  < 3 • <19, 

where  7ft)  ~ lOOt  hours. 

Pfoblaa  lr  Compute  E„(r)  and  E.^(t)  for  9 • 0.  9. . a.  9_,  ard  “O  „ 

Solution  From  table  i wa  can  raad  Ee(r)  for  this  case,  sinca  k “ eo/01^  • 3 

and  a - - o05.  For  9 - 0,  0^  (-  2500),  s (-  4115),  9q  7500),  oO 

respect  ■sly  ws  haw  BQ(r)  * 3,  6.1,  7.2,  2,9,  0.  Eg(t)  la  found  moat 

vun'Ily  : the  replacement  caaa  for  all  values  J 9{j*oo)  by  using  (12), 

EgCt)  - EgCr'".  Thus  wa  gat  for  9 • 0,  9^,  a,  9Q  respectively,  EQ(t)  - 0, 

1 57" a ^ -20,  Far  9 « <J«  , tha  expected  waiting  tin*  to  roach  a decision 

-lOOt^/3753  . 

is  give,  by  t*,  » ahera  a - ~ * This  gives  - E ^(t)  - 110, 

--9q  log  B 

Remark:  'rore  gana rally,  in  tama  of  B,  n,  9 , and  k,  t • „ 

n(k  - 1) 

Thi3  raa  -3  that  if  no  items  fail  by  time  t ^ stop  experimentation  at  t rjJ 
with  »c<  aptance  of 

Problem  Assume  that  we  are  besting  the  hypothesis  in  Problem  1.  A 
sample  ’ sics  100  is  placed  ot,  test.  Items  which  fail  are  replaced  by 
new  it«j  > drawn  from  the  same  lat.  The  experiment  ia  started  at  ti-Jr  t - 0* 
Th-  fir  failure  occurs  at  x-  *■  20~1  hours,  the  second  failure*  ?t 

•A.  * 

Xg  - 1C  5 hours,  th*  third  failure  at  Xj  * 121.7  ho’ir?. 


tha  fourth  failure 
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at  * io'/o4  tours,  the  filth  failure  at  x^  - 175.2  hours.  (Ail  lii-^s  are 
measured  from  t *»  0. ) 

(a)  Verify  that  rs>  decision  has  bean  reacted  by  tias  x^. 

(b)  Wa  keep  waiting  for  the  sixth  failure  and  note  that  it  has  net  yet 
occurred  at  287,5  hours  (tiae  measured  from  t - 0).  Verify  that  we 
can  atop  experimentation  at  tine  t * 267*5  with  the  acceptance  of  H .. 

Solution:  It  can  be  readily  verified  that  in  this  case  (5)  becosse 

-100  ♦ 37. 5r  < t < 100  ♦ 37. 5r.  This  region  ie  drawn  in  Figure  1.  The 
life  test  data  are  plotted  on  the  figure  b/  sorug  vertically  so  long  as 
v»  are  waiting  for  the  next  failure  to  occur  end  moving  horizontally  by 
one  unit  (in  r)  at  each  failure  time.  Clearly  the  path  crosses  into  tha 
region  :f  acceptance,  when  r - 5,  at  tine  t ” 100  * (37*5)5  “ 28?., 5*  Sine* 
thu  sixth  failure  has  not  yst  occurred  we  can  stop  experimentation  at 
t • 287,5  with  the  acceptance  of  H#. 

Rsisarkt  ku  a matter  of  fact  we  happen  to  know  in  this  ease  that  the  sixth 

failure  occurs  at  x*  - 346.7  hours.  Thus  we  s^ved  346.7  - 237*5  • 59=2  hours 

■*>* 

by  observing  the  life  trst  continuously  in  time. 

Fr5bl«a  j:  Ths  first  seven  failure  tis*s  in  a -ample  of  100  (with  replace- 

msnt)  are  x,  - 19.  3;  Xg  ••  45.8,  « 49.9,  ■ 96.7,  - 115.2, 

■ 12?s?,  *-  - 121.2.  Verify  that  if  the  hypotheses  be3;«g  tested  are 
thos*  in  Problem  1,  then  is  at  tias  x^  *»  131.2  hours. 

Solution ; See  Figure  20 

Remark:  Not?  that  while  the  decision  in  Problem  3 is  mads*  between  x^  and 
the  decision  in  Problem  4 i»  made  at  the  failure  time  Xy.  Also  iter#  is  an 
excess  over  the  boundary. 
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Pr»falcn  5t  Find  a truncated  ( non- sequent  ial ) refinement  procedure  for 
testing  the  hypothesis  in  Prnblsre  3L  Uee  constant  sample  size  n - ICO* 
Solution:  From  results  in  ^6_7,  it  can  be  verified  that  the  truncated 

replsucsnscit  procedure  nesting  tbs  requirmsnts  i.3  aa  foils'**: 

If  i*0? o $J  - 407. 5,  truncate  the  experiment  at  407.$  with 

acceptance  of  H#.  If  min  jx^,  407.51  • x,q,  truncate  sacDerlmentation  at 
x^  with  acceptance  of  The  0.0.  curves  of  this  test  procedure  and  the 
one  in  Preb'  -x*  1 are  eseertialiy  the  cameo 

Prvhl—  6:  Ccespute  Eq(t)  and  E^(t)  for  the  plan  in  PrwbXea:  5 .ter 

9 ■ 0,  Sj.,  a-  9^,  96  , 

Selutfon:  From  results  in  /~6_7,  Eft(r)  - 10,  9.93,  8.75,  5.39t  0, 

E0(t)  - | ®e(**)o  For  9-0,  9,,  a,  9#,  oo  respectively,  Eg(t)  - 0,  248, 
360,  404.5,  407.5  respectively. 

Remark:  In  Figures  3 and  4 we  coopare  the  EQ(r)  and  E0(t)  curves  for 
Problems  2 and  6,  This  will  giv*  some  idea  of  the  savings  in  the  expected 
number  of  failures  and  tin*  to  reach  a decision. 

Problss  7:  Find  in  Froblf*»  1 if  a « & - .01. 

-0„  log  B 

Solution:  t__  = — »■*  * 230.  This  is  about  twica  th-  value  of  t _ 

n(k  - 1) 


when  a - P - .05. 


VfrlUBE 


in  Hour® 


V- 


r » nuaber  of  failures. 

1 

I Fi^urs  Is  Graphical  tr«aira*ni  of  dais  in  rr-^ilaa  3, 


Ki 2 lira  Zi  Gi4.jiu.C4Ll  traataant  of  data  In  Probl*  A 
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